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Abstract 

The stabihty of the optical pulse of the Crab pulsar is analyzed based on the 1 fis resolution observations with the Russian 6-meter 
O^and William Hershel telescopes equipped with different photon-counting detectors. The search for the variations of the pulse shape 
^ along with its arrival time stability is performed. Upper limits on the possible short time scale free precession of the pulsar are 

placed. The evidence of pulse time of arrival (TOA) variations on 1.5-2 hours time scale is presented, along with evidence of small 

light curve (shape and separation of main and secondary peaks) changes between data sets, on time scale of years. Also, the fine 

structure of the main pulse is studied. 
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1. Introduction 

■ Over the last 30 years the Crab pulsar has been ex- 
tensively studied. The reasons for it are clear - it is the 
brightest pulsar seen in optics, it is nearby and young. 
However, the most popular groups of contemporary the- 
ories of the Crab high-e nergy emission, the "polar cap" 
llDatighhertv et al . 19961 ) and "outer gap" ( Cheng et al . 
ones, can't explain the whole set of observational 

data. 

One of the main properties of the Crab emission is the 
relatively high stability of its optical pulse shape despite 
the secular de crease of the luminosity, related to the spin 
rate decrease ( PaciniL 1971 : Nasuti et a 

At the same time the pulsars in general and the Crab it- 
self are unstable. The instabilities manifest themselves as 
glitches, likely related to the changes of the neutron star 
crust, timing noise, powered by the collective processes in 
the superfluid internal parts of it, magnctospheric instabil- 
ities, results of the wisps around the pulsar, precession, etc. 
All these factors may influence the optical pulse structure 
and change it on various time scales, both in periodic and 
stochastic way. 

However, it has been found early that the variations of 
the Crab optical light curve, in contrast with the ra dio ones , 
are g overned by the Poissonian statistics (Kris tian et all 
19701) . A number of observations show the absence of non- 



stationary effects in the structure, intensity and the dura- 
tion of the Crab optical pulses, and the restrictions on the 
regular and stochastic fi ne structure of its pulse on the time 
scales from 3/iS to 500/is feeski n et al . 1983 : Percival et al 



1993 ) , the fluctuations of the pulse intensity (jKristian et all . 
1970|). 



Along with the increase of the observational time spans 
and the accuracy of measurements, small changes of the 
optical pulse intensity, syn chronous with the giant radio 
pulses, have been detected (jShearer et all . Sooi). Also, the 
evidence for the short time scale precession of t he pulsar has 
been found by studying its optical light curve ( Cadez et aj 
200 ih . 

All this raises the importance of monitoring the Crab 
optical emission with high time resolution. 

The article is organized as follows. In Section 2 we briefly 
describe the observation process and instruments used, in 
Section 3 the method of phase stability analysis is described 
and used to study the Dec 1999 and Jan 2007 data sets, 
in Section 4 the light curves of different data sets are com- 
pared, in Section 5 the possible fine structure of the main 
pulse peak is discussed, and Section 6 gives the conclusions. 

2. Observations and data reduction 

We analyzed the sample of observational data obtained 
by our group over the time span of 12 years on different tele- 



Table 1 

Log of observations 



Date 


Telescope 


Instrument 


Duration 
seconds 


Spectral range 

A 


Dec 7, 1994 


BTA, Russia 


Four-color photometer 


2400 


U -1- B -1- V -1- R 






with photomultipliers 






Dec 2, 1999 WHT, Canary 


Avalanche photo-diode 


6600 


4000-7500 




Islands 








Jan 9, 2000 


BTA, Russia 


Panoramic photometer 


7900 


U -1- B -1- V -1- R 






with position-sensitive detector 






Nov 15, 2003 


BTA, Russia 


Avalanche photo-diode 


1800 


4000-7500 


Jan 25, 2007 


BTA, Russia 


Panoramic photometer 


10500 


B -1- R 


Jan 26, 2007 


BTA, Russia with position-sensitive detector 


6500 


B -f R 



scopes. The details of observations are summarized in Table 
[H The equipment used were a four-color standard photome- 
ter with diaphragms based on photomultiplie rs, a fast pho- 
tometer with avalanche photo-diodes (Shearer et al , l2003l) 
and a pa noramic photometer based on position-sensitiv e 
detector (IPebur et all l2003l : IPlokhotnichenk o et all . l2003l ). 
All devices are photon counters which record the photon 
time of arrivals with accuracy better that at least 1/iS, and 
the final observational data are t he lists of these t imes. 
In case of panoramic photometer ( Debur et al . 20031 ) only 
photons arrived in 3" aperture around the pulsar location 
has been considered. 

The photon lists of all observational sets have been pro- 
cessed in the same way by using the same software to ex- 
clude the systematic differences due to data analysis in- 
consistencies. Photon arrival times of 1999, 2003 and 2007 
years sets have been collected with absolute time scale cal- 
ibration by means of GPS receivers. So, they have been 
corrected to the barycenter of the Solar System using the 
adapted version of axBary code by Arnold Rots. The ac- 
curacy of this code has b een te sted with detailed examples 
provided by Tvne et a 1 (|2003h and is found to be better 
than at least 2^s. 

The barycentered photon lists then have been folded in- 
dependently u sing both Jodrell-Bank radio ephemerides 
( Jordanl[200i) and o ur own fast-folding base d method of 



timing model fitting ( Plokhotnichenko . 1993h . 



The declared accuracy of Jodrell-Bank ephemerides fre- 
quency and derivatives provide the folding precisi on of at 
least 1 IIS, and the base epoch - of at least 5 /iS (| Jordan . 
l2006l) . so we decided to fold the light curves with 5000 bin 
(6.6 s) resolution. 

The observations of 1994 and 2000 have been per- 
formed without absolute time scale calibration, using non- 
stabilized frequency generators, so these data can't be con- 
verted to the barycenter correctly. So, we divided the whole 
data set into several pieces short enough to be well fitted 
with a 3rd order timing model (up to second frequency 



derivative), performed an independent timing model fit for 
each, folded and combined it together, compensating the 
phase shift between separate pieces. The accuracy of such 
procedure is proved to provide similar time resolution, so 
we use the same number of bins in its analysis. 

The observations of 1994 and 2000 have been performed 
in standard Johnson-Cousins U, B, V and R photometric 
bands, 2007 - in B and R bands, while 1999 and 2003 - 
without fil t ers, b ut with the same detector, described by 
Ryan et al ( 20061 ) (avalanche photodiode with broad spec- 
tral sensitivity in the 4000-7500A peaked at - 7000A). 

3. Study of phase stability 

We performed the search for timing model residuals 
( "phase shifts" ) using two longest continuous data sets of 
1999 and 2007 years. 

3.1. Computation of timing residuals 

The photon list of 1999 data set has been divided into 
segments of 100 pulsar periods (with duration of ap- 
proximately 3.3 s), which resulted in 1677 segments with 
mean number of photons of 10546.5 each. They then have 
been folded separately using the same Jodrell-Bank ra- 
dio ephemerides, and each fold have been cross-correlated 
with the template, which has been built by folding the 
entire data set. The estimation of the segment phase shift 
in respect to template is then derived by approximating 
the peak of cross-correlation function (in a phase window 
of 0.02 period width) with the Gaussian, and analytically 
computing its maximum position. The steps of the proce- 
dure are illustrated in Fig. [TJ The resulting phase shifts 
are plotted in upper panel of Fig. [3l 

The formal accuracy of the maximum of Gaussian ap- 
proximation of cross-correlation function estimation is 
much better than the spreading of real data values (basi- 
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Fig. 1. Method of estimation of phase shift of the sample folded 
light curve in respect to the template one. Upper panel - sample 
light curve of 1999 data set with 5000 bins resolution and 100 pulsar 
periods long (approximately 3.3 sec). Superimposed is the template 
profile, derived by folding the whole data set. Middle panel - cross- 
correlation of the sample and the template light curves, with Gaus- 
sian fit superimposed. Lower panel — cross-correlation residuals after 
subtraction of Gaussian fit. The real accuracy of phase shift estima- 
tion is much worse than the one of determining the cross-correlation 
peak position due to influence of original light curve errors and its 
systematic deviation from the Gaussian approximation. However, it 
may be shown that it does not lead to the statistical biasing of the 
estimate. 




Number of samples 

Fig. 2. Dependence of the estimation of phase shift mean value of 
simulated data on the set length. The simulated data have the same 
parameters (template profile, number of lightcurve bins and mean 
number of photons per sample) as the 1999 set. No significant biasing 
is seen, and the deviation from zero decreases with the increase of set 
length. It proves that the estimation of phase shifts is statistically 
unbiased and reliable. 

cally, the error bars are hidden inside the dots in Fig [3]). 
This is partly due to neglecting of original light curve er- 
rors while computing the cross-correlation. Also, the lower 
panel of Figure [1] demonstrates that the Gaussian estima- 
tion for the cross-correlation function is not perfect, as it 
shows systematic deviations from it. To ensure that these 
facts are not spoiling the results and to test the statistical 
quality of the method used we performed the numeri- 



Fig. 3. Upper panel — phase residuals of the Crab pulsar after applying 
the third-order timing model (up to second frequency derivative). It 
corresponds to the Gaussian noise with 5.5 ■ 10~* cycles rms. Lower 
panel — results of the similar analysis of the simulated data with the 
same mean parameters. The statistical properties are roughly the 
same as for the real data set. 



Optical data, 1999 - 
Radio data, 1982-1989 (Scott et ai 2003) ■ 




Frequency, 1/day 

Fig. 4. Power density spectrum of the timing residuals of the 1999 
data set, computed according to method descri bed in | Sc ott et a] 
l|2003l ) . Also, the spectrum of radio residuals from lScott et all l|2003h 
is shown along with its power-law fit with a slope a = —3.09. 

cal simulation by generating the set of sample Poisson- 
distributed light curves based on the average profile of 
the 1999 data set with the same mean number of photons 
per segment, and processing them in the same manner 
as the real data. The computed phase shifts of simulated 
data are shown in the lower panel of Fig. [3l The RMS 
of the simulated and observed data are roughly the same 
(cr Ri 5.5 • 10^"*), which is consistent with the statistical na- 
ture of phase shift values scatter. Also, we performed the 
test whether the estimated phase shifts are unbiased and 
statistically reliable by studying the behaviour of the mean 
value of simulated data phase shift and its RMS in depen- 
dence of number of segments. The results of this simulation 
are shown in Fig. O It may be easily seen that the mean 
value of the estimation along with its RMS, both converge 
to zero with the increase of the number of segments, which 
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Fig. 5. Power density spectrum of the timing residuals of the 1999 
data set, computed without any time domain window. The horizontal 
line shows the mean noise level. The deviation from it of the first 
bin value has the significance level SL = 2.4 • 10""* and suggests the 
presence of variations on time scale comparable with data set length. 

proves the statistical reliability of the method used. Also, it 
permits to increase the determination of systematic phase 
shifts accuracy by averaging it over long phase segments. 



3.2. Fourier analysis of phase shifts 

In order to search for periodic components of the phase 
shift we performed the Fourier anal ysis of the data se t ac- 
cording to the method described in lScott et all (j2003r ). i.e. 
computed the power-density spectrum P using the time- 
domain Hann window Wi oc sin^ {-ni/N) to suppress the 
power leakage (which lowers the spectral resolution approx- 
imately by two) . 

The resulting power density spectrum is shown in 
Fig. Ulin comparison with results of radio data analysis of 
Scott et al ( 2003 ). The accuracy of our data is not suffi- 
cient to reach the level of extrapolated power-law timing 
noise seen in the radio band. 

However, it is possible to derive upper limits for the si- 
nusoidal variable components of timing noise on 3.3 s - 50 
minutes time scale. Indeed, for the purely white noise pro- 
cess the val ue 2P/ (P) is dist ributed as with 2 degrees 
of freedom ( Leahv et all . ll983h and has mean and standard 
deviations of 2. The probability Q for this quantity to ex- 
ceed some threshold value 2Pq / (P) by chance is 



Xo 



2^ 



(1) 



For the significance level SL = 0.01 (which corresponds 
to the 99% confidence probability) the threshold value Pq, 
determined by solving 



2Po 

(P) 



Fig. 6. Phase shifts of the 1999 year data set (upper panel, each 
point corresponds to 3.3 seconds of data) and its averaging in 10 
time segments 562 seconds long (lower panel). The vertical error 
bars correspond to the standard deviations of values in segments. 
Also, the approximation of phase shifts with 3-rd order polynomial. 
Characteristic time scale of variations is 1.5-2 hours with amplitude 
~ (1.5 ± 0.5) ■ 10-* (4±1.5 ^is). 

is Pq — 4.6 (P) , where (P) is the mean "noise" level of power 
density spectrum. By combining it with the spectral ampli- 
tude Pa — ^ (where Sv = 1/T is the spectral resolution) 
of sinusoidal signal, we have 



2Ai/ 



< 4.6(P) 



and for the upper limit for amplitude 



A < v/9.2(P)Ai. 



(3) 



(4) 



(2) 



For our data, the mean "noise" level of power density is 
(P) = 2.2 • 10"" cycles^-day and Aj/ = 15.3 day-\ which 
leads to the limit for amplitude of periodic component of 
timing residuals A < 5.6 • 10~^ cycles (1.8 /is) in this fre- 
quency range. 

3.3. Secular behaviour of phase shifts 

To test the phase stability of the Crab light curve on time 
scales comparable with the total length of observations we 
first computed the power density spectrum of phase shifts 
without any time domain window. The result, shown in Fig- 
ure[5l demonstrates significant deviation (with significance 
level SL = 2.4-10"'') of the first bin from noise mean value, 
which suggests the presence of variations on time scale of 
the data set length. 

To check these variations in time domain we divided the 
data set into 10 equal time segments 562 seconds long, and 
computed the mean and variance of phase shifts in it. The 
results (shown in Figure [6]) indeed show the presence of 
significant variations on 1.5 hours time scale with (3 ± 1) • 
10^** (9±3 ^s) amplitude. 

We performed similar analysis of the data of Jan 25- 
26, 2007. The signal to noise ratio in this set is smaller, 
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Fig. 7. Phase shifts of the Jan 25, 2007 data set (upper panel, each 
point corresponds to 33 seconds of data) and its averaging in 10 
time segments 1050 seconds long (lower panel). The vertical error 
bars correspond to the standard deviations of values in segments. 
Also, the approximation of phase shifts with 7-rd order polynomial 
is shown. Characteristic time scale of variations is 1-2 hours with 
amplitude ~ (3 ± 1) ■ 10"*' (9±3 ^ts). 
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Fig. 8. Phase shifts of the Jan 26, 2007 data set (upper panel, each 
point corresponds to 33 seconds of data) and its averaging in 5 time 
segments 1300 seconds long (lower panel). The vertical error bars 
correspond to the standard deviations of values in segments. Also, the 
approximation of phase shifts with 3-rd order polynomial is shown. 

SO we computed the phase shifts using 1000 period long 
segments of hght curve (roughly 33 seconds). The results 
(see Figures [7] and [8]) also show variations on 1.5-2 hours 
time scale with similar amplitude. 

4. Pulse shape 

We performed the comparison of pulse profile shapes 
of data sets of 1994, 1999, 2000, 2003 and 2007 years in 
a wa.y similar to the one used in ( Jones. Smith fc Nelsonl . 
19801) . As there a re evidences that t h e Crab p rofile depends 
on w ave length 
(|200nh : iBeskin. 



Eikenberrv et al.' (1996'); 'Golden et al. 
Komarova fc Plokhotnichenka (2000): 



z_ 
-500 r- 



Fig. 9. Crab profile change in B and R filters between 1994 and 2000. 
Intensity is in arbitrary units - all light curves are normalized to the 
Nov 2003 data set one. R-band data demonstrates the combination 
of pulse phase shift (~ 30-50 fis) and small change of peaks right 
wings. 



Romani et all (|2001l ). we compared the data acquired in 
the same wave bands, i.e. 1994 vs 2000 and 2007, and 1999 
vs 2003 data sets. 

We compensated the phase shifts of 1999 and 2007 data 
described above by approximating it with high-order poly- 
nomials and adding it to the timing models. Also, to com- 
pensate possible phase shifts between light curves of dif- 
ferent data sets (due to, for example, systematic errors in 
radio ephemerides base epochs) we determined the phase 
shift between them by the method described in Section [3] 
and re-folded them with base epoch shifted according to 
it. This procedure has been performed iteratively until the 
phase shift became smaller than at least lO^'*, i.e. less than 
half of a bin size used. 

Then we normalized each light curve i/i to the same tem- 
plate profile yo^i (we used the one of Nov 2003 data set as 
it has the largest number of photons) by means of linear 
transformation y'^ = ayi + b with parameters a and b max- 
imizing the likelihood 



L'^Y. (y.lnA, -ln(y,!)-A,) , 



(5) 



i=0 



of yt to be the instance of Poissonian distribution with 
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Fig. 10. Crab profile change in B and R filters between 2000 and 
2007 data sets. Both bands demonstrate peaks (primarily — first one) 
shape change without significant phase shift. 

1,, b 

Then we rebinned them in blocks with roughly equal 
number of photons and plotted the difference between 
them. The results are presented in Figures [9lfT2l 

To check whether it may be due to uncompensated phase 
shift we simulated shifted light curves, computed their dif- 
ferences and plotted them in Figure[6l The simulated effect 
has different shape and different ratio of positive and neg- 
ative residuals, which argues for the reality of the detected 
variation. 

The effects seen in FigsCTT^mav be interpreted as a some 
combination of systematic "phase shift" , variation of main 
peak shape and change of the distance between primary and 
secondary peaks. Unfortunately, its exact nature can't be 
revealed by the methods used - it is impossible to correctly 
define the "phase shift" of two profiles with different shape, 
and the procedure of light curves phasing devours some 
part of the shape change effect. However, it may only lower 
the significance of detected residuals - so the presence of 
the effect itself is undoubtful. 

5. Pulse fine structure 

For the data set of Nov 2003, which has the largest num- 
ber of photons collected, we performed the search for the 



Fig. 11. Crab profile change in B and R filters between 1994 and 
2007 data sets. Change of front wings of peaks is combined with the 
decrease of peak separation. 
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Fig. 12. Crab profile change between 1999 and 2003 data sets. The 
data are unfiltered, but acquired on the same photometer. Phase 
shift is combined with main peak right wing steepening. 

fine structure of the peak of main pulse. Wc approximated 
its shape by means of low-pass Fourier filter with character- 
istic frequency of 0.05 bin~^ , which effectively smoothes the 
light curve with sinc-like window with ~ 40 bins FWHM. 
The peak data, fit and residuals are shown in Figure [T31 
We do not detect any significant spike-like fine structure 
on the level of 0.5% (f a) with 6.6/xs time resolution. 
However, there is a single "absorption-like" feature with 
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Fig. 13. Simulated residuals of the identical light curves with sys- 
tematic shift of 1-5 bins (6.6 - 33 fis). The shape of residuals is not 
like the ones discovered in comparison of different data sets. 
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Fig. 14. Approximation of the peak of the Crab pulse profile of Nov 
2003 data set by means of Fourier low-pass filter with characteristic 
time scale of 20 bins (132 /is). The diagonal cross marks the feature 
with A = 3.77(T amplitude. The probability of its appearance on the 
peak (above the half-amplitude of the main pulse, 227 bins of 5000) 
by chance is 1.9%. 

an amplitude of S.77a near the maximum. The probabihty 
for such a feature to appear by chance on the peak (on 
the segment with intensity above the half-amphtude which 
contains 227 bins) is 1.9%. 

The hght curves of other data sets, as well as their sum, 
however, do not show any statistically significant devia- 
tions. 



6. Discussion 

Optical emission of Cra b pulsar has been d e tected 
soon after its radio o n e (IStaelin fc Reifensteinl . 1X968 : 
Cocke. Disnev fc Tavlon . 119691 ). an d since then i t has 
been observed a number of times dKristian et al 1970 [ 



Cocke & FergusonVl974';'Peterson et al.','l978';'Percival et al.', 
1993;iBeskin. Komarova fc Plokhotnichcnka,200Q;, Golden et 



2000f) ). Being the brightest (-16™) among 5 known 
optical pulsars, it demonstrates relatively high stabil- 
ity of light curve shape agai nst a b a ckgro u nd predicted 
secul ar luminosity changes (|Pacinil I1971I: iNasuti et all . 
19961) . At the same time, it exhibits the timing noise in 



pulse time of arrivals in radio, optical and x-ray spectral 
bands on time scales from se v eral days to tens o f years 



(ICorded. Il980l: iBovnton et all Il972t iKuiper et all 12003 : 



Rots. Jahoda fc Lvnd . 2004 ). Also, the variation of pri 



mary and secondary peaks intensity ratio in gamma rays 
has given the hint of 13-year periodicity (jNolan et al.l . 
1993t Ulmer et al.l . Il994 ) . Similar behaviour has b een ob- 



serve d in radio on months and years time scales 
1988t iLvne. Pritchard fc Smithl . Il988t IScott et all . 




On the other hand, the search for time of arrival residuals 
on short (seconds to hours) time scale has not been practi- 
cal ly performed. Th e one except ion seems to be the result 
of ICadez fc Galici3 (fl99Q' ) and ICadez et all (|200ll ). who 
detected the 60-seconds periodicity of pulsar intensity. 

We analyzed the data of several sets of optical observa- 
tions with high temporal resolution of the Crab pulsar per- 
formed by our group over the last 12 years. 

No evidence for periodic short time scale variations of 
pulse time o f arriv als (like 60-sec free precession claimed by 
Cadez et"^ ( 2001 )) is detected on the 3.3 s - 50 minutes 



time scale on Dec 2, 1999. The upper limit for their am- 
plitude is ^ < 5.6 • 10~^ cycles (1.8 fis) (significance level 
0.01). Note that no periodic features has also been detected 
in the Crab light curve on these frequencies in previous 
work of iGolden etal\ (|2000l ). 

Also, no signature of extended timing noise spectral fea- 
tures (like power-law one observed on lower frequencies) is 
seen on this time scale. 

However, the data of Dec 2, 1999 and Jan 25-26, 2007 
sets both show significant phase variations on 1.5 - 2 hours 
time scale with (2 — 5) • 10"*^ cycles (6 - 16 /xs) amplitude. 
This effect is most likely not truly periodic. Moreover, it 
is difficult to explain it as a precession of a rotating rigid 
body, as it requir es too large d i fferen ce of neutron star axes 
(AR/R - 10-5) (|Akgun et a]| . l2006h . 

Possible manifestations of noise processes, related to 
either superconducting vortices inside the neutron star, or 
magnctospheric effects, on the short time scales of minutes 
to hours has yet to be analyzed. These effects has been 
usually involved in explanatio n of timing noise observed on 



time scales of days to years (ICordes fc Greensteinl . 11981 
lAlpar et aill986l : lcheng et all Il987ai bh. However, the am- 



plitude of the effect we discovered significantly exceeds 
the power- l aw ex trapolation of timing noise spectrum 
(jScott et all . l2003l ) (see Figgj). 

Only one observed effect is currently kn own to occur on 
simil ar time scale - the giant radio pulses ( Lundgren et aj 
Il995h . which have inverse power-law intenstity distribu- 
tion and randomly appear in all phases of light curve oc- 
cupied by "normal" radio emission except for precursor 
('jes sner et at [20051 ). Moreover, it has been recently sug- 
a^^sted that all radio emission except for the precursor con- 



sists of giant pulses only (jPopov et all . l2006l ) . Their origin 
is most likely due to changes of either coherence conditions 
or electron density in the magnetosphere. In the latter case, 
it may influence the optical emission region too. The slight 
correlation between giant radio pulse s and increase of op - 
tical emission has been discovered in ( Shearer et a 
As giant pulses appear randomly in phase, they may lead 
to changes of optical pulse shape, and so mimic the time of 
arrival variations. 

Also, it may in principle be attributed to polar cap 
current-pattern drifti ng, which may occur on a very broad 
range of time scales (jRuderman fc Gil . 2006f) . 

The non-detection of this effect in ongoing radio observa- 
ti ons may be attributed to i t s lowe r resolution (according 
to Lvne. Pritchard fc SmithI ( 1993h . the accuracy of pulse 
time of arrival determination in Jodrell Bank observations 
is 20 /iS for 10 min integration time). 

We discovered the variation of pulse shape between 
different sets of our observations, i.e. on time scale 
of several years. It presents and has similar proper- 
ties in all studied spectral bands, and cannot be at- 
tributed to well-known effec t of s h ape dependence on 
wavelength (lEikenberrv et all. I l996': 'Golde n et all. I2OOOI: 
Beskin. Komarova fc Plokhotnichenko. 2000; iRomani et al . 
200 Due to limitations of data analysis methods used it 
is impossible now to specify the exact nature of the varia- 
tion - it may only be empirically described as a combina- 
tion of systematic phase shift, main and secondary peaks 
shape change and variation of peak separation. Also, it is 
not clear whether the variation periodic, systematic or ir- 
regular. Ho wever, the effect is similar to th e one marginally 
detected in ( Jones. Smith fc Nelson . 1980[ ) on time scale of 
7 years. 

There are several possible physical mechanisms able 
to produce such pulse shape variations on time scale of 
years. First is the susp ected precession of Crab on ^ 568 
days ( Scott et "al . 2003h . Indeed, at least one other pulsar 
- PSR B1828-11 - exhibits the precession accompanied 
by th e changes of a radio pu l se pro file on a similar time 
scale ( Stairs. Lvne fc Shemar . 2000[ ). For Crab, however, 
due to difference in rotational frequencies, such preces- 
sion period implies much smaller wobble angle, and so - 
smaller pulse profile variations. Also, polar cap current- 
pattern drifting may mimic the precession and result in 
the same phase shift and p r ofile change effects on years 
time scale. ( Ruderman fc Gil . l2006h . 

Another possibi lity is th e incomplete post-glitch relax- 
ation (Demiaiiski fc Proszvhs ki. 1983; L vne. Pritchard fc Smit^ 
1993t IWong. Backer fc Lvnd . l200li) . as all our observa- 
tions have been separated by glitches of different power 
( Lvne et"al . [2003) . In Crab, it manifests itself as a per- 
sistent change of frequency derivative, and may be at- 
tributed to small changes of the angle between the mag- 
netic dipole and the ro tation axis (Link fc Epstein, 1997; 
Allen fc Horvathlll997tlLink et a]| . ll998h . which inevitably 
leads to pulse profile variation. Also, pulse profile variations 
in hard energy band are often observed in anomalous x- 



ray p ulsars ( Kaspi et al. . 2003t Morii. Kawai fc Shibazaki , 
20041) ■ however, it is still not clear whether they result 
directly from glitches. 

We do not detect any spike-like fine structure of the main 
pulse maximum on the level of 0.5% (1 a) with 6.6fis time 
resolution. 

All the proposed explanations of discovered variations of 
pulse shape and time of arrival are qualitative only and are 
in no sense complete. The observations have to be contin- 
ued, and the theoretical analysis still has to be performed. 
We hope the study of such variations can help to elaborate 
the theory of pulsar emission. 
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